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We present numerical studies of non-linear conduction in graphene nanoribbons when a bias 
potential is applied between the source and drain electrodes. We find that the conductance quan- 
tization plateaus show asymmetry between the electron and hole branches if the potential in the 
ribbon equals the source or drain electrode potential and strong electron (hole) scattering occurs. 
The scattering may be at the ends of a uniform ballistic ribbon connecting wider regions of graphene 
or may be due to defects in the ribbon. We argue that, in ribbons with strong defect scattering, the 
ribbon potential is pinned to that of the drain (source) for electron (hole) transport. In this case 
symmetry between electron and hole transport is restored and our calculations explain the upward 
shift of the conductance plateaus with increasing bias that was observed experimentally by Lin et 
al. [Phys. Rev. B 78, 161409(R) (2008)]. 

PACS numbers: 72.10.Fk,73.23.Ad,72.80.Vp,73.50.Fq 



I. INTRODUCTION 



Conductance quantization in graphene nanoribbons 
was recently observed experimentally by Lin et al. Un- 
like in the ballistic quantum wires in which conductance 
quantization had been observed previously, 2 the quan- 
tized conductance steps observed by Lin et al. 1 were 
two orders of magnitude lower than the conductance 
quantum 2e 2 /h. The observation 1 of such low quan- 
tized conductance values was surprising. However, it 
was subsequently shown theoretically 3,4 that it can arise 
from enhanced electron backscattering near subband 
edges of disordered graphene ribbons with carbon atom 
vacancies 3 or chemisorbed atoms or molecules. 4 More re- 
cently, conductance quantization has also been reported 
in low conductance ribbons that were fabricated in a 
different way, 5 and ballistic conductance quantization 
in graphene nanoconstrictions has also been observed. 6 
With increasing bias voltage applied between source and 
drain electrodes the quantized conductance plateaus were 
observed by Lin et al. to shift upward for both electron 
and hole conduction. It was suggested 1 that this effect 
was due to the increase with increasing bias of the num- 
ber of subbands that are present in the energy window 
between the source and drain electrochemical potentials 
and therefore contribute to transport. However, whether 
this mechanism can by itself explain the experimentally 
observed non-linear transport effect 1 has not been inves- 
tigated theoretically and is still an open question. In this 
paper we investigate non-linear conduction in graphene 
ribbons with the help quantum transport simulations. 
We find that strong electron scattering by defects and 
a specific type of potential profile in the ribbon under 
bias are both necessary for the the nonlinear transport 
phenomenon observed by Lin et al. to occur. 

We describe the models of the graphene ribbons and 
potential profiles that we study in Section II. We present 
our results for ballistic ribbons having uniform widths 
and having wide-narrow-widc geometries in Section III. 
Our results for ribbons with defects are presented in Sec- 



tion IV A and the physics of the potential profiles in such 
ribbons is discussed in Section IV B. Our conclusions are 
summarized in Section V. 



II. MODEL 

We describe graphene ribbons by the standard tight- 
binding Hamiltonian on a honeycomb lattice, 
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where ej = Vu is the on-site energy that accounts for 
bias potential Vf, at site i; tij = t = 2.7 eV is the matrix 
element between nearest-neighbor atoms. This Hamil- 
tonian is known to describe the it band dispersion of 
graphene well at low energies. 7 Spin and electron interac- 
tion effects are disregarded in present study. The effect of 
disorder is represented by bulk vacancies that are intro- 
duced by randomly removing carbon atoms and setting 
appropriate hopping elements ty to zero. It is quantified 
by the probability of the carbon atoms being removed 
p, which is normalized relative to the whole sample. We 
refer the reader to our previous papers, Ref. 3,4, and 
other studies 8-10 for the effects disorder of other types 
on conduction in graphene ribbons. 

Electron transport out-of-cquilibrium is modeled by 
application of a source-drain voltage V sc i that biases the 
electrochemical potentials of the source and drain leads 
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spectively, where E-p is the Fermi energy. Note that 
throughout this paper our definitions of the source and 
drain are such that electrons flow from the source to the 
drain. The potential profiles within experimentally re- 
alized graphene ribbons are not well understood since 
in addition to the source-drain and gate voltages they 
depend on the presence of charged adsorbates on the 
graphene, charged defects in the substrate and other un- 
charactcrizcd materials properties. We will therefore con- 
sider three basic model potential profiles that might be 



2 




Energy (t) 

FIG. 1: (color online) Conductance as a function of the Fermi 
energy for ideal uniform (a) and QPC (b) graphene ribbons. 
The gray solid lines show zero bias conductance V B d = 0; the 
green dotted lines correspond to linear potential ramp along 
the structure; the solid blue and dashed red lines correspond 
to the entire potential drop V a d occurring at the source or 
drain end, respectively, of the ribbon in (a) and constriction 
in (b). The gray dotted line in (b) shows the conductance of 
the uniform ribbon for comparison; it is the same as the gray 
solid line in (a). The width and length are 10 nm and 100 
nm for the ribbon in (a) and constriction in (b); the leads in 
(b) are twice as wide. The edges are armchair except for the 
wide-narrow boundaries in (b) that are shown in the inset. 
The arrow in (b) marks the Fermi energy for the DOS plots 
in Fig. 2. t = 2.7 eV. 



realized in experimental devices and compare the results 
of our quantum transport calculations for these model 
profiles with the experimental results of Lin et al. 1 The 
profiles that we will consider are: a linear ramp of T4, 
and constant Vb along the channel with a potential drop 
either at the source or drain electrode. We will show 
here that only potential profiles of the latter type are 
consistent with the nonlinear transport data of Lin et 
al. 1 and will discuss the physical origin of such potential 
profiles in Section IV B. For a non-zero V s d the electric 
current is calculated as 11 I = ^ dE T(E, V s( i) with 
T(E, V s d) being the transmission coefficient. Here and 
below we assume zero temperature. In order to calculate 
T{E 1 V s d) we employ the non-equilibrium Green's func- 
tion method, 12,13 which is also used for calculation of the 
density of states (DOS) in the ribbon. The conductance 



of the ribbon is given by the Landauer formula 13 
J Op 1 

G= V^ = TTT dET(E 7 V sd ). (2) 

Vsd n Vsd J^d 

For graphene ribbons with vacancies we also present the 
averaged conductance (G) that was obtained by averag- 
ing over 10 samples with different realizations of disorder. 

III. BALLISTIC RIBBONS 

Let us consider first the non-linear conductances of 
ideal graphene ribbons of finite length L with the source- 
drain potential applied between its two ends. The latter 
are attached to semi-infinite leads that serve as electron 
reservoirs and have constant electrochemical potentials. 
The bias potential is allowed to drop inside the region 
of length L. Because Vj, is not known apriori we ex- 
amine three generic situations: i) linear potential ramp 
along the ribbon, ii) sharp potential energy drop near the 
source, iii) sharp potential energy drop near drain elec- 
trode. For simplicity the potential changes in two latter 
cases are assumed to occur over a distance less than the 
lattice constant. The results presented here are for the 
armchair edge orientation of the ribbons. Note that sim- 
ilar results were obtained for zigzag edges as well. 

A. Uniform Ribbons 

Fig. 1(a) shows the conductance of an ideal ribbon of 
uniform width for different V s d and bias potential profiles 
as a function of the Fermi energy. The ribbon length 
and width are chosen to be L = 100 nm and W = 10 
nm, respectively This corresponds to 235 unit cells and 
N = 80 carbon atoms in transverse direction of the 
armchair edge-oriented ribbon. 14 Note that such ribbon 
sizes are readily accessible using present day fabrication 
techniques. 1,15 The ideal leads that carry electrons to and 
from the ribbon are represented by semi-infinite uniform 
graphene ribbons of the same width. The main effect 
of the source-drain bias applied to the ideal ribbon that 
is seen in Fig. 1(a) is gradual broadening (that is sym- 
metrical about the Dirac point E = 0) of the transitions 
between the quantized conductance plateaus. This is a 
result of the finite energy window eV s d that contributes 
to the current in Eq. (2). Within that window the trans- 
mission probabilities T due to states belonging to differ- 
ent subbands are averaged in Eq. (2) resulting in less 
steep transition regions between quantization plateaus. 
Note that a similar effect occurs in quantum point con- 
tacts (QPCs) in the two-dimensional electron gas. 16 The 
conductance for a given V s d is determined by availability 
of particular initial states for tunneling and their trans- 
mission through the channel. Only states in an energy 
window of width |e|V" sc i = /j, 3 — /i d participate in transport. 
In the case of ballistic transmission, the integral in (2) 
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becomes (n-l)(|e| V sd -e') + ne' if e' = fJ, s -E n < |e| V sd , 
where E n is the energy of the bottom of the highest occu- 
pied subband n in the source. Thus, as [i s changes rela- 
tive to E n , which might be realized by changing the Fermi 
energy, the conductance (2) varies gradually between the 
n — 1 and n quantization plateaus making the transition 
region eV s d wide, see Fig. 1(a). When e' > |e|T4d the 
integral in (2) will be n and the conductance reaches a 
plateau. Another feature of electron transport in uniform 
ideal graphene ribbons is the insensitivity to the details 
of the potential profile that is evident in Fig. 1(a): Very 
little reduction of the conductance occurs if the smooth 
linear potential profile is replaced by a sharp potential 
drop in the ribbon; the small change might be attributed 
to increased electron scattering at the potential discon- 
tinuity. However, the quantized conductance values in 
Fig. 1(a) clearly do not increase as the bias voltage in- 
creases and more subbands contribute to transport, in 
marked contrast to the behavior that was observed ex- 
perimentally by Lin et al. 1 Thus the fact that the number 
of subbands that contribute to transport increases with 
increasing source drain bias is not by itself sufficient to 
explain the nonlinear transport phenomenon that Lin et 
al. 1 observed. 



B. Wide-Narrow- Wide Ribbon Geometries 

We now consider a QPC-like structure that is more 
relevant to experimental devices, 1,15 where the graphene 
ribbon connected wider areas of graphene. Fig. 1(b) 
shows the conductance of such a graphene ribbon QPC 
with the wider regions twice as wide as the constriction, 
TV = 158 and N = 80 carbon atoms in cross-section, re- 
spectively. The constriction and wider regions both have 
armchair edges that arc connected by zigzag edges at the 
wide-narrow boundaries, see the inset in Fig. 1(b). Both 
of these widths support a propagating state at zero en- 
ergy, which means both the wide and narrow armchair 
regions are metallic. 14 The dimensions of the constric- 
tion are similar to those of the ribbon studied in Fig. 
1(a). Despite this similarity the conductances in Fig. 
1(a) and (b) are quite different. This is due to strong 
electron scattering at the interfaces connecting the nar- 
row and wide regions of the QPC ribbon. Such strong 
scattering was discussed in Refs. 17,19 and attributed 
to extra electron degrees of freedom on the hexagonal 
graphene lattice. Note that electron transport here dif- 
fers substantially from that in conventional QPC's in the 
two-dimensional electron gas in GaAs heterostructures: 
No adiabatic transmission occurs in the graphene ribbon 
QPC's even if the interfaces between the narrow and wide 
regions are made atomically smooth. 1 The zero bias con- 
ductance of the QPC ribbon shows strong oscillations 
whose envelope resembles the step-like behavior of the 
uniform ribbon; compare the solid and dotted lines in 
Fig. 1(b). The oscillations result from Fabry- Perot reso- 
nant scattering at the entrance and exit of the ribbon. 18 




FIG. 2: (color online) Representative LDOS and DOS for the 
QPC ribbons for the Fermi energy Et — 0.07i and source- 
drain bias V s d = 0.025t/|e|. DOS was calculated by integra- 
tion across the ribbon. The constriction potential equals the 
source (a) or drain (b) potential. The small outsets on the 
right show the transmission coefficient. 



They are quickly smeared by bias if eV s d exceeds the 
oscillation period and quantization steps become clearly 
discernible in that case. As the bias increases further the 
quantization steps are smoothed in a way similar to the 
uniform ribbon in Fig. 1(a). However, two important 
features appear in addition. First, the conductance is 
suppressed more strongly, see the green dotted lines in 
Fig. 1(b) for \e\ V sd = 0.025t, 0.05i, 0.075t where the lin- 
ear bias ramp was assumed. In addition to the widening 
of the transition regions between plateaus, the plateaus 
themselves are lowered due to poor transmission through 
the QPC constriction. Second, the non-linear conduc- 
tance is asymmetric about the Dirac point if the model 
electron potential energy in the constriction Vb is such 
that the entire potential drop V s d occurs at one end of 
the constriction. For example, if it occurs at the end con- 
nected to the source (so that Vb = — \ |e| V s d is "pinned" 
to the drain) then the conductance curve shifts to neg- 
ative Fermi energes; see the blue curve in Fig. 1(b) for 
|e| V s d = 0.025<. The size of the shift is proportional to 
Vsd- The opposite shift occurs for Vb pinned to the source. 
These asymmetries, however, become quickly damped for 
large V s d due to conductance suppression (not shown). 

The conductance asymmetry in Fig. 1(b) that is trig- 
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gered by electron scattering at the ends of the QPC con- 
striction can be understood by considering the DOS plots 
in Fig. 2. The figure presents cases of source (a) and 
drain (b) "pinned" constriction potentials for a represen- 
tative Fermi energy Ep = 0.074. Conduction through the 
QPC structure is constrained by the constriction, not by 
the electrodes where there are more transport states at 
every energy. If the constriction potential equals that 
in the source, Fig. 2(a), electrons from the second and 
third subband in the leads have no corresponding avail- 
able states in the constriction and transmission in the 
energy range eV s d surrounding Ep is damped. The to- 
tal transmission is close to 1 due to the first subband 
state that still propagates freely in the constriction at 
E = 0.074. Inspection of the LDOS in the upper panel 
of Fig. 2(a) reveals low amplitudes inside the constric- 
tion. If Vb = — |e| Vsd/2 as in Fig. 2(b) the second and 
third subbands in the constriction become available for 
electron transmission to the drain electrode. The trans- 
mission is, however, not perfect due to state mismatch at 
the interfaces. 17,19 The DOS inside the constriction shows 
multiple resonant peaks due to the Fabry-Perot interfer- 
ence. Each peak is due to a resonant state available for 
propagation with enhanced transmission T at the corre- 
sponding energy. The LDOS at the Fermi energy for the 
drain-pinned bias potential is much higher than for the 
source-pinned potential; compare the upper plots in Fig. 
2(a) and (b). The larger LDOS is a manifestation of the 
higher probability to find an electron in the constriction. 
Note that by far the strongest electron localization occurs 
at the corners of the structure regardless of the potential 
distribution. If the Fermi energy changes sign and the 
charge carriers become holes the availability of transport 
states inside the constriction reverses. This results in the 
conductance curve being shifted in energy by an amount 
proportional to V s d relative to its position at zero bias. 
The conductance shows no asymmetry about the Dirac 
point for the linear ramp potential because of approxi- 
mately similar availability of transport states inside the 
constriction for the Fermi energies of opposite sign. 

In Fig. 1(b), for the QPC models with the poten- 
tial drop at one end of the constriction, the conduc- 
tance increases with increasing source drain bias in a 
limited range of low bias values for some values of the 
Fermi energy. However, these conductance increases oc- 
cur in the transition regions between the quantized con- 
ductance plateaus. The quantized conductance values on 
the plateaus themselves do not increase with increasing 
bias although the Fermi energy (or gate voltage) ranges 
in which the plateaus occur shift as the bias increases. 
This behavior is again qualitatively different than that 
observed by Lin et al. 1 who reported the quantized con- 
ductance values increasing with increasing bias without 
the ranges of gate voltage in which the quantized plateaus 
occurred changing significantly. We conclude that the ob- 
served nonlinear transport phenomenon 1 is unlikely to be 
a property of ballistic graphene ribbon QPC's. We will 
show next that electron scattering due to disorder can 
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FIG. 3: (color online) The averaged conductance vs. Fermi 
energy for 30 nm wide and 500 nm long graphene ribbon with 
p = 10 -4 bulk vacancies. The dashed and solid lines corre- 
spond to \e\ V 3 d = 0.0124 and 0.0244, respectively. Note that 
valley degenerate subbands are separated by ~ 0.0244 for 30 
nm wide ribbon; the subband positions are marked by verti- 
cal dashed yellow lines. The blue and red lines show results 
for source and drain pinned constriction potentials. The gray 
filled area denote zero bias conductance. The branches of the 
conductance curves that resemble the experimental non-linear 
conductance data of Lin et al. 1 are outlined by the thick dark 
gray dashed and solid lines. Averaging is performed over 10 
samples with different realizations of the bulk vacancies. 



give rise to this effect. 



IV. UNIFORM RIBBONS WITH DEFECTS 
A. Results 

We calculated the averaged conductances for uniform 
ribbons with bulk vacancy disorder. Representative re- 
sults are shown in Fig. 3. The ribbon width was chosen 
to be W = 30 nm, the same as that of the experimental 
samples in Ref. 1. The presence of the bulk vacancies 
leads to enhanced electron backscattering near subband 
edges that in turn suppresses the conduction equally for 
all subbands resulting at low bias in conductance steps 
of equal height 3 similar to those observed by Lin et al. 1 
Edge disorder as well as long-range potentials are proba- 
bly present in the experimental samples 1 but have been 
shown 3 not to be of prime importance for the observed 
conductance quantization and we, therefore, concentrate 
here on scattering by bulk vacancies. It is worth not- 
ing that adsorbed H atoms cause electron scattering and 
conductance quantization that are very similar to those 
due to bulk vacancies. 4 

Figure 3 shows the calculated non-linear conductances 
for the ribbon with bulk vacancies when the bias po- 
tential in the ribbon is pinned to the source or drain 
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potential, i.e. V b = +\e\V sd /2 or V b = -\e\V sd /2. 
The plots are for two values of |e| V s d = 0.012i and 
0.024i that differ by half of a subband energy separa- 
tion. When |e| V s d increases from 0.012i and 0.024i the 
conductance plateaus shift upward (with some smooth- 
ing) by half of the conductance quantization step for elec- 
tron transport if V& = — |e| V s d/2 and for hole transport 
if Vb = + |e| Vsd/2. The corresponding branches of the 
conductance curves are outlined by the thick dark gray 
dashed and solid lines in Fig. 3 and resemble the experi- 
mental results in Fig. 4 of Ref. 1. That is, in the present 
model, the conductance plateaus move to higher conduc- 
tance values with increasing bias as in the experiment 
of Lin et al. 1 if the potential in the ribbon is pinned to 
that in the drain for electron transport and to that in the 
source for hole transport. By contrast, if a linear ramp 
potential in the ribbon is assumed, the quantized conduc- 
tance plateaus are found not to shift to higher conduc- 
tance values with increasing source-drain bias. Instead 
the transitions between them broaden until the plateaus 
disappear. 



ribbon the potential is expected to be pinned at the drain 
potential for Fermi levels above the Dirac point (electron 
carriers) and at the source potential for Fermi levels be- 
low the Dirac point (hole carriers), in agreement with 
the model potential profiles that we have shown above to 
yield behavior of the nonlinear quantized conductances 
resembling that observed experimentally by Lin et al. 1 

The nonlinear conductance in the experiment 1 exhib- 
ited plateaus at V s d = 100 mV with conductances higher 
than those for V sd = 50 mV by roughly one unit of the 
quantized conductance but with similar plateau widths. 
Note that 50 mV is energy separation between valley- 
degenerate subbands. However, our results show the 
quantization steps becoming smoother as the bias in- 
creases. This discrepancy might be attributed to the ef- 
fect of electron-electron interactions. At the same time 
the experimental paper 1 presented just a single measure- 
ment for a single sample. Figure 3 shows, by contrast, 
a statistically averaged result. Therefore, sample-specific 
conductance quantization should not be ruled out when 
experimental data in Ref. 1 is analyzed. 



B. Discussion 

The model ribbon potential profiles that we have con- 
sidered above have been phenomenologically motivated. 
However, we shall now argue that it is physically reason- 
able that Vb is approximately uniform throughout most 
of the ribbon and that it changes from being pinned to 
the drain to being pinned to the source when the Fermi 
energy crosses from above to below the Dirac point: Cal- 
culations of the potential profiles in idealized models of 
short (a few nanometers long) gated graphene ribbons 
have indicated that due to screening most of the potential 
drop in the ribbon that is due to the applied source-drain 
bias occurs near the contacts. 20 The characteristic length 
scale for such potential variations has been estimated 21 
to be A = ^(e GNR /e x)tGNRtox where e GNR and e x 
are the graphene ribbon and gate oxide dielectric con- 
stants and tcNR and tpx are the graphene ribbon and 
gate oxide thicknesses. For the ribbon of Lin et al. 1 this 
yields A ~ 10 nm which is much smaller than the ribbon 
length of 1.7 /mi in their non- linear conductance experi- 
ment. This suggests that our assumption that Vb is con- 
stant throughout most of the ribbon is justified. Wc also 
note that since the measured conductances of the ribbons 
of Lin et al. 1 were much smaller that 2e 2 /h most of the 
charge carriers entering the ribbon were reflected back 
into the electrode from which they came without pene- 
trating deeply into the ribbon. Thus most of the charging 
of the ribbon by the carriers and hence most of the vari- 
ation of the electrostatic potential in the ribbon should 
have occured near the end of the ribbon where the carri- 
ers enter, i.e., near the source electrode for electrons and 
the drain for holes. Consequently throughout most of the 



V. CONCLUSIONS 

In conclusion, we have studied conductance quanti- 
zation in graphene ribbons in the non-linear transport 
regime. We found that with increasing source-drain 
bias the transitions between the quantized conductance 
plateaus broaden. We also found that the conductance 
plateaus show significant asymmetry between the elec- 
tron and hole branches if the potential in the channel is 
pinned either to the source or to the drain electrode po- 
tential and strong electron (hole) scattering occurs. This 
scattering may occur at the ends of a uniform ballistic 
ribbon that connects wider regions of graphene or may 
be due to defects in the interior of a ribbon. We argue 
that in a ribbon with strong defect scattering the poten- 
tial in the ribbon should be pinned to the drain potential 
for electron transport and to the source potential for hole 
transport. In that case the electron- hole symmetry of the 
conductance is recovered and our results explain the up- 
ward shift of the conductance plateaus with increasing 
bias between source and drain observed experimentally 
by Lin et al. 1 for both electron and hole transport. Fur- 
ther theoretical study that accounts for electron-electron 
interactions is required and more experimental data are 
needed to better understand non-linear electron trans- 
port in the graphene ribbons. 
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